Microscopic Study of the triple-a Reaction 
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We present a microscopic dynamical study of the reactions involving three "^He clusters. We show 
that the much discussed triple-a linear chain configuration of ^^C is formed with a certain lifetime 
and subsequently decays into a triangular configuration of ^^C and then to a configuration near the 
ground state. Time-dependent Hartree-Fock (TDHF) theory coupled with a density constraint is 
used to study the properties of these configurations. We find a sequence of dynamical transitions 
analogous to the suggested astrophysical mechanism for the formation of ^^C nuclei. 

PACS numbers: 21.60.-n,21.60. Jz,21.30.Fc,21.60.Cs,27.20.+n,27.30.+tl 



Reactions involving three "^He nuclei carry both astro- 
physical and nuclear structure significance. Astrophys- 
ically, the triple-a reaction is believed to be one of the 
dominant helium burning reactions leading to the syn- 
thesis of most of the ^^C found in the universe P, 0- 
In the first stage of this process two helium nuclei com- 
bine to form an unstable ^Be nucleus, which has an ex- 
tremely short lifetime (~ 10~^^s) and decays into two 
q:'s. However, for densities attained during the helium 
burning process there will be a small abundance of *Be 
present to allow the combination with another "^He nu- 
cleus to form the well known 3] resonant state of ^^C 
at 7.654 MeV, which is at an energy close to the com- 
bined energies of the reacting ®Be and ^He nuclei, thus 
facilitating this reaction 2]. This is the reason why de- 
spite its non-existence under natural conditions ®Be plays 
a crucial role in nuclear astrophysics. The formed reso- 
nant state of ^^C is believed to decay to the 2+ state of 
^^C at 4.43 MeV and subsequently to the ground state 
via 7 radiation. Recently, there has also been interest 
in non-resonant reaction of three alpha's as an alternate 
mechanism at lower temperatures [4|. 

The exotic structure of the state of ^^C and its 
possible association with the mysterious state of ^^O 
at 6.06 MeV has also been of great interest in nuclear 
structure models of light nuclei. In 1950's it has been 
conjectured that the state might be linear chain con- 
figuration of 3 alpha's [5[ , however later it has been inter- 
preted to be gas-like, consisting of a mixture of various 
configurations of three alpha's including the linear chain 
configuration Therefore, search for the linear chain 
configurations of alpha's has been performed in different 
places on the nuclear chart. Suggestions that the life- 
time of the linear chain configuration may be extended 
by the addition of extra neutrons have led to the discus- 
sion of cluster structures in neutron-rich Be 0-11] and C 
isotopes [13" 12 1 • The search has also been extended to 



other light N = Z nuclei |13l - ll6l |. However, until now 



there has been no theoretical discussion on the dynamics 
of the linear chain state; e.g. how the chain state of three 



alpha's plays a role in the early stage of the synthesis of 
^^C or what are the characteristic modes of vibration of 
the three alpha chain. Another problem is that, until now 
most of the theoretical analyses for the cluster structures 
have been performed using effective interactions, which 
are determined such as to reproduce the binding energies 
and scattering phase shifts of the clusters. It is highly de- 
sirable to study the presence of cluster configurations in a 
more general manner without the a priori initialization in 
terms of clusters. Also, dynamical oscillation of the linear 
chain configuration is a completely new approach to the 
problem. In our previous study, cluster configurations 
of neutron-rich C isotopes have been obtained as local 
minima of static Hartree-Fock (HP) solutions with ax- 
ial symmetry, and TDHF calculations were performed to 
examine their stability against quadrupole and octupole 
distortions [TtI . [18| . The linear chain configurations were 
found to be relatively stable with respect to the breathing 
(quadrupole) distortions while they were unstable with 
respect to bending (octupole) perturbation. 

It is generally acknowledged that the TDHF theory 
provides a useful foundation for a fully microscopic many- 
body theory of low-energy heavy- ion reactions [l9| . Ear- 
lier TDHF calculations of collisions between light nuclei 
involving cluster-like configurations have been made for 
the study of nuclear molecular resonances [13] . However, 
due to the lack of computational resources these calcula- 
tions suffered from numerical imprecision as well as un- 
physical symmetry assumptions, such as collisions being 
restricted to axial symmetry. Current TDHF calculations 
are performed with high numerical precision and with no 
symmetry assumptions as well as using modern Skyrme 
forces. Recently, we have shown that when TDHF is com- 
bined with the density-constraint method (21j dynamical 
potentials and ion-ion interaction barriers can be accu- 
rately reproduced 22l424l| . Thus an alternative approach 
would be to investigate the formation and stability of the 
linear chain configuration using the fully microscopic and 
dynamical TDHF theory. 

In this work we study TDHF collisions which repro- 
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duce the linear chain configuration and subsequently de- 
cay to lower-energy configurations of the system. To our 
knowledge such mode changes have never been observed 
in TDHF calculations previously and appear to simulate 
the suggested astrophysical mechanism for the formation 
of ^^C nuclei. 
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FIG. 1: (Color online) Selected density profiles from TDHF 
time-evolution of the ''He-l-^Be head-on collision for initial Be 
orientation angle /3 = 0° using the SLy4 interaction. The 
initial energy is -Ec.m. = 2 MeV. 

In our numerical calculations we have chosen a Carte- 
sian box which is 40 fm along the collision axis and 24 fm 
in the other two directions. Calculations are done in 3-D 
geometry and using the full Skyrme force (SLy4) [2^ as 
described in Ref. [26[. We have tested these results for 
few cases by using a much larger numerical box and by 
changing the time-step and found no appreciable differ- 
ence. Similarly, using different parametrizations of the 
modern Skyrme force all show the same phenomenon. 
We have chosen to study two different collisions lead- 
ing to the linear chain configuration, first the ''He-l-^Be 
system and then the triple collision of three '*He nuclei, 
which may be astrophysically much less probable. The 
Hartree-Fock (HE) state for the ^Be nucleus is axially 
symmetric. In Fig. [T]we show three snapshots from the 
long time evolution of the ''He-l-^Be collision at zero im- 
pact parameter and Ec.m. = 2 MeV. The top panel of 
Fig. [T] shows the linear chain configuration about which 
the system oscillates for times less than 2, 500 fm/c. In 
particular, it is remarkable that the moving clusters do 
not equilibrate while moving inside the linear chain state 
but retain their 2p — 2n character, where one observes a 
complex quasiperiodic motion with little damping up to 
this time. Around 2500 fm/c the system starts to bend 
and acquires a somewhat triangular shape as shown in 
the middle panel of Fig. [TJ The system still retains its 



cluster character with the center cluster moving off the 
reaction plane cut shown in the figure, but can be clearly 
observed in volumetric three-dimensional movies of the 
collision process. The bending motion, where the cen- 
ter cluster oscillates somewhat perpendicular to the left 
and right clusters continues for approximately 1000 fm/c, 
with very little damping. Finally, at even longer times 
the system relaxes into a relatively more compact shape 
(bottom panel of Fig. [T]). Such mode changes, where the 
dynamical energy in the longitudinal direction is con- 
verted to a transverse mode, while the system retains its 
cluster structure have never been seen in previous TDHF 
calculations albeit this would not have been possible in 
calculations imposing axial symmetry. Even in three- 
dimensional calculations, for an exactly central collision, 
the axial symmetry would be preserved under ideal theo- 
retical and numerical conditions. Therefore the meaning 
of head-on collision {b = 0) should be interpreted to have 
a small dispersion around this value, which facilitates the 
mode change even for exactly central collisions. 
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FIG. 2: (Color online) Single-particle parities of the neutron 
states during the collision of the ''He-l-^Be system as a func- 
tion of time at -Ec.m. ~ 2 MeV. 

We were also successful, for the first time, in creat- 
ing a static Hartree-Fock linear chain state orthogonal to 
the ground state in a fully three-dimensional geometry. 
This was achieved by initializing one of the single-particle 
states to be in the s — d shell with positive parity rather 
than in the p state with negative parity. This results in 
a linear chain state similar to the one shown in the top 
panel of Fig. [TJ The fact that this state has two positive 
parity and one negative parity state proves that it is ex- 
actly orthogonal to the ^^C ground state which has one 
positive parity and two negative parity states. A similar 
dependence on parity was also studied in cluster model 
calculations [27| . In order to relate this observation to the 
dynamical mode changes discussed above we have used 
the density constraint method to calculate the potential 
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energy and the single-particle parities during the TDHF 
time-evolution. In Fig. [2] we show the neutron single- 
particle parities as a function of collision time. What 
is striking is that the combined system initially has the 
same parity signature as the static linear chain state but 
at the time of bending one of the positive parity states 
starts to decay towards negative parity and this decay 
continues as the system becomes closer to the parity sig- 
natures of the ground state. Oscillations in the numer- 
ically calculated parities stem from the fact that these 
are done during the dynamical evolution of the system. 
The proton single-particle parities are almost exactly the 
same as for the neutrons as anticipated. 
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FIG. 3: (Color online) Potential energy curves for the collision 
of the "^He-l-^Be system as a function of R for three initial 
alignments of the Be nucleus and at -Ec.m. = 2 MeV. 

In order to relate the observed mode changes more 
closely to the intrinsic energy of the system we have also 
calculated the potential energy of the system as a func- 
tion of the ion-ion separation distance R 22| . For the cal- 
culation of R we have used the hybrid method described 
in 2J] , which relates R to the quadrupole moment of the 



system thus making it possible to have a consistent def- 
inition of R for large overlaps. The calculated potential 
energy curves are shown in Fig. |3] as a function of three 
alignments of the ^Be nucleus with respect to the colli- 
sion axis labeled as angle /3, and for the entire duration of 
the collision process. Since for the real system the angu- 
lar momentum of ®Be is projected all possible alignments 
of the *Be nucleus needs to be considered. For all of the 
alignments the combined system climbs up a shallow po- 
tential barrier height of approximately 1.24 — 1.44 MeV, 
the lowest barrier being that of the (5 — Q° potential, 
thus making this alignment most probable under thresh- 
old conditions. For the potential energy curve showing 
the head-on collision (/3 — 0°, black curve) we observe 
that the system initially relaxes to a relatively shallow 
metastable minimum and oscillates about this minimum 



until approximately T = 2500 fm/c at which points it 
slips down the curve towards the second configuration as 
indicated by three alphas in a triangular configuration. 
After spending some time in this configuration the sys- 
tem further slips down to even lower energy and more 
compact configuration. The potential energy curves cor- 
responding to ^Bc initial alignment angles of /? = 45° and 
(3 — 90° (red and blue curves, respectively) undergo a dif- 
ferent behavior, bypassing the linear chain minimum but 
directly going to the triangular and subsequently to the 
compact configuration, the perpendicular energy collision 
attaining the most compact and lowest energy configu- 
ration. It is interesting to note that all of the potential 
energy curves spend some time in the triangular config- 
uration. 

As an alternate collision leading to the same config- 
uration we have also studied the collision of three ''^He 
nuclei, one at rest at the origin of the collision axis and 
the other two on each side boosted towards the center 
with 1 MeV energy. In Fig. U] we contrast the time de- 
pendence of the potential energies for the two different 
collisions. We observe that the three ''He collision process 
spends considerably longer time (about 6000 fm/c) un- 
dergoing quasiperiodic oscillations with very little damp- 
ing in the linear chain configuration before switching to 
bending and compact modes. 
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FIG. 4: (Color online) Time development of the potential 
energy for the head-on collision of the *He-|-®Be and the 
''He+''He-f-''He systems for Ec.m. = 2 MeV. 

In order to gauge the stability of the linear chain con- 
figuration we have made systematic studies as a function 
of impact parameter and center-of-mass energy, as well 
as a study using heavier Be isotopes to determine the de- 
pendence on neutron number. In Fig. [5] we show the de- 
pendence of the linear chain survival time on the impact 
parameter for the ''He-l-^Be system at Ec.m. — 2 MeV 
and (3 = 0°. We observe that as the impact parameter in- 
creases the survival time rapidly decreases. This decrease 
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naturally happens slower (faster) for lower (higher) en- 
ergies. We have also done a similar study for the time 
spent in the linear chain configuration as a function of the 
ccntcr-of-mass energy for the "'He+^^Be system for /3 = 0° 
alignment. We decreased the energy in steps of 0.1 MeV 
to find the lowest energy for which we form the linear 
chain configuration (at lower energies the nuclei rebound 
due to Coulomb repulsion). At this energy of 1.3 MeV 
the lifetime of the linear chain configuration increases to 
about 2875 fm/c. As the energy is increased the lifetime 
decreases gradually. In order to study the dependence of 
the linear chain state on the neutron number of the Be 
isotopes we have repeated all of the above calculations 
using a ^Be nucleus instead. The calculations were done 
by using all the time-odd terms in the Skyrme interac- 
tion appropriate for an odd-A nucleus. While we do find 
an analogous behavior in this study, the lifetime of the 
linear chain state is approximately 30% less than that of 
the corresponding ^Be system. 
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FIG. 5: Time spent in the linear chain configuration as a 
function of the impact parameter b for the ''He+^Be system 
at Ec.m. = 2 MeV and /3 = 0° alignment. 

The analysis of the long-time motion of collective nu- 
clear phenomena in terms of their quasiperiodic behavior 
is related to the long-time behavior of nonintegrable me- 
chanical systems and addresses questions as to the nature 
of energy dissipation and equilibration of energy. In this 
work we have performed microscopic dynamical calcu- 
lations of nuclear collisions to study the formation of a 
metastable linear chain state of ^^C. The time evolution 
as calculated using the TDHF equations shows a char- 
acteristic quasiperiodic exchange of alpha-like clusters in 
the density function corresponding to a quasiperiodic mo- 
tion along a static Hartree-Fock potential energy surface, 
which is studied using the density constraint procedure. 
We have shown that the calculations lead to the forma- 
tion of a metastable linear chain state of three alpha-like 
clusters which subsequently decays to a lower energy tri- 



angular alpha-like configuration before acquiring a more 
compact shape. This is the first observation of such mode 
changes in TDHF calculations and the results seem to be 
analogous to the suggested astrophysical mechanism for 
the formation of ^^C nuclei. 
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